1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.03
## LDevsum 1.01 1.05
## dh0 1.38 2.18
## dh1 1.72 3.73
## dh2 1.47 2.65
## dl0 1.09 1.34
## dl1 1.19 1.68
## dl2 1.03 1.03
## dl3 1.19 1.66
##
## Multivariate psrf
##
## 1.38
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1227.65640 | 22087.6761 |
| DIC3 | 1160.16296 | 22151.3761 |
| PWAIC | 44.21523 | 242.1802 |
| WAIC | 1187.89048 | 22173.1181 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.04
## LDevsum 1.08 1.33
## dh0 1.41 2.23
## dh1 2.51 5.00
## dh2 1.57 3.22
## dl0 3.28 6.91
## dl1 1.07 1.08
## dl2 3.35 7.44
## dl3 1.76 3.22
##
## Multivariate psrf
##
## 3.05
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1381.02334 | 22212.1649 |
| DIC3 | 1246.20600 | 22234.6104 |
| PWAIC | 90.86349 | 297.1184 |
| WAIC | 1315.67881 | 22266.9946 |